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ON  THE  STABILITY  OP  THE  ROTATIhiG  MOVEMENTS  OP  A  SOLID  BODY  THE  CAVITY 
OP  WHICH  IS  PILLED  WITH  AN  IDEAL  LIQUID 

by 

N,  G,  Chetayev 

The  problem  of  the  aocviracy  of  the  flight  of  a  liquid-filled  missile  along 
its  trajectory  is  a  hard  one  if  the  liquid  is  vlsoous.  If  we  consider  the  li¬ 
quid  in  the  missile  as  ideal  and  incompressibley  then  we  can  find  the  ccrreot 
solutions  of  the  problems  in  the  stability  of  rotating  movements  of  this  solid 
body,  provided  that  the  interior  is  completely  filled  with  liquid  without  any 
bubbles.  In  my  opinion  this  solution  does  not  seem  unimportant,  since  starting 
from  it  we  can  consider  as  suffioient  the  surplus  of  stability  against  the  nega¬ 
tive  influences  of  visoosity. 

1,  General  Suggestions.  Among  the  many  possible  displacements  of  a  mis¬ 
sile  filled  with  an  ideal  Inoompresslble  liquid  free  from  bubbles,  are  the  fol¬ 
lowings  rotation  around  any  straight  line,  and  the  forward  displaoements  like 
that  of  a  solid  body.  Hence  the  movement  of  this  missile  is  governed  by  the 
theorem  on  moment  in  a  system  with  its  origin  in  center  of  gravity  of  the  mis¬ 
sile  and  the  liquid  0,  and  with  the  axes  ?  k  ^  parallel  to  rigid  axes. 

f  I  ^ 

These  conditions  allow  us  to  exaimine  in  the  problem  of  the  stability  of 
the  rotating  motions  of  a  missile  whose  interior  is  filled  with  an  ideal  liquid, 
only  the  relative  motions,  considering  the  oenter  of  gravity  0  as  being  fixed. 

The  axes  x,  y,  z  are  Introduced  to  facilitate  the  calculations  let  z  be 
the  axis  of  the  elliptical  Inertial  rotation  of  the  missile  (without  liquid) 
plotted  at  point  0,  perhaps  not  coinciding  with  the  oenter  of  gravity  of  the 
solid  shell;  and  let  the  axes  x  and  y  lie  in  the  plane  perpendicular  to  the  z 
axis,  as  may  be  more  convenient  for  us,  but  in  such  a  manner  that  the  moments 
of  inertia  of  the  solid  shell  with  reference  to  the  three  axes  will  always  be 
_ the.oonstant  and  main  axes  of  the  ellipsoid  of  the  inertia  of  the  shell  plotted 

at  point  0. 
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It  is  assumed  that  the  motion  of  the  ideal  Incompressible  liquid  filling 
the  cavity  of  the  missile  Is  determined  by  the  Instantaneous  velocities  of  the 


solid  shell.  Let p  be  the  constant  density  of  the  liquid  ;  n,  outside  normal 
to  the  surfaoe  S  of  the  cavity  occupied  by  the  liquid;  and  direc¬ 

tional  cosines  of  the  normals  on  the  axes  x,  y,  z.  Let  p,  q,  r,  stand  for  In¬ 
stantaneous  angular  velocities  of  the  rotation  of  the  missile,  around  the  axes 
X,  y,  z. 

At  the  beginning  of  motion  the  liquid  was  stationary  In  the  missile,  liy 
the  Lagrange  theory  the  liquid  will  move  with  the  potential  of  the  velooities^^ 
(x,  y,  z,  t). 

The  function  ^  is  determined  under  the  condition  that  the  normal  compon¬ 
ent  of  the  velocity  of  the  liquid  coincides  with  a  point  on  the  shell  in  such 


a  manner  that 


=  (<;z  —  ry)  a  -  j  (rx  —  pz)  ft  -(  ( PV  -  i 


Let  us  assume 


and  let  the  relationship  (l)  be  satisfied  independently  from  the  terms  p,  q,  r. 


This  gives 


%=y7--*P. 


Henoe  ,  does  not  depend  on  t,  but  depends  on  the  cavity  which  Is  filled 
with  an  Incompressible  ideal  liquid. 

The  kinetic  energy  cf -the  rotating  motion  of  the  missile's  shell  Is 

2T' =  Ay C'r^ 

where  A',  B',  C  are  the  constant  moments  of  Inertia  of  the  missile's  solid 
shell  together  with  the  dividing  walls  of  the  cavity  with  respect  to  axes  i,  y,  z. 
Kinetic  energy  of  relative  motions  of  the  liquid,  in  the  cavity 
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2^-=fni[(5)'+(2r)‘+©‘]‘''=  s 

“1  =  /’.  '"t  =  7<  = 

II  ,t.  I— „ fff fiii'i ^  j- IrfT 

{'!'<•  ^/1  ^JJJ  \_  dz  dx  '  ay  dy  '  a  di  } 


2,  Cavity  in  the  Form  of  a  Circular  Cylinder.  Let  us  examine  the  cavity 
in  the  form  of  a  circular  cylinder  with  radius  a,  ajcis  z,  and  height  2h.  Since, 
in  respect  to  any  straight  line,  passing  through  the  point  0  and  perpendicular 
to  the  axis  z,  the  moment  of  Inertia  of  the  hard  shell  has  onl^'’  one  value,  it 
is  easier  to  determine  the  position  of  the  missile  by  using  the  coordinate  sys- 
tem£)^  j",  connected  with  the  missile,  (Figure  l).  (see  page  1?) 

The  system  of  rectangular  axes  (£hf)  when  rotated  around  the  axis  n  at  an 
angle  <X  jvill  change  into  the  system  (  but  the  system  )  when  rota¬ 

ted  around  the  axis  x  at  an  angle changes  into  the  system  (xyz).  The  angle  of 
rotation  of  the  missile  around  its  axis  z  in  the  system  (xyz)  we  will  denote 
by  OJ  . 

Prom  the  diagram  we  nwst  oonolude  that  the  angles,  are  holonomous 

ooordinates  of  the  missile,  and  hence  the  differential  equation  of  the  rotating 
motions  of  the  missile  with  these  variables  will  resemble  Lagrange  equation. 

Prom  the  diagram  we  conolude  that  the  projeotion  of  the  absolute  instant¬ 
aneous  velooity  of  the  rotating  solid  shell  along  the  axes  x,  y,  z  are 

P  =  ^>  y  =  a'co3p,  r  =  —  a'sinp-f***’ 

Where  a',  p',  «'  are  derivatives  of  time  t  corresponding  to  the  angles  «.  P,  w. 

Total  kinetic  energy  of  the  missile  together  with  its  liquid  in  the  inter¬ 
ior  will  be  T  ■  T'  +  T<*. 

Similarly  to  the  classic  example  of  the  stability  of  the  rotating  motions 
of  a  oonventional  shell  with  solid  equipment  on  a  flat  trajectory,  let  us  con¬ 
sider  only  the  resisting  oouple  with  moment  ,  proportional  to  the  sine  of 

the  angle y  between  the  z  axis  of  the  missile  and  the  velooity  of  the  center  of 
gravity  in  the  system  0, 
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Assuming  that  the  axis  C  is  directed  along  velocity  of  the  oenter  of  grav¬ 
ity  in  system  0,  then  we  shall  have  '  cos  7  =  cos  a  cos  p. 

The  possible  work  of  the  resisting  couple  is  simply 

|i  sin  7J7  =  ^( — cos  7)  =  i{ — |iCos*cos  P)  = 

=  sin  a  cos  p8a  + 1*  cos  a  sin  p8p  si  ^,8a  +  ^^8p  + 


From  this  the  generalized  forces  Qa  are  determined  by  the  formulas 

^•=  (1  sin  a  cos  p,  Qp=(icosasinP,  ^^  =  0 

In  order  to  set  up  the  differential  equations  of  rotation  motions  of  the  missile 
filled  lyith  liquid,  one  must  calculate  T*.  The  problem  of  determining  the  func¬ 


tion  'p  was  solved  by  N,  E,  Zhukobskiy,  The  solution  is  given  below  with  insig¬ 
nificant  amplification. 


The  boundary  condition  of  function  on  the  sitrfaoe  S  of  the  hollow  gives 
d^Jdn  =  0,  since  on  the  lateral  surface  of  the  cylinder^ *-^^and  on  the  bases  of 
the  cylindera  =  0,  p  =  0.  Here  by  Neumann's  principle  we  have  i|(,=  const,  and  conse¬ 


quently 

lij.  'l'sl=0.  IV  'l's)=0»  fV  '{'»)  =  0 


Prom  the  ajcial  symmetry  of  the  cavity  it  oan  be  concluded,  that  it  is  suf¬ 
ficient  to  determine  only  one  of  the  functions  ;  <{«,,  <{«,.  .  Let  us  determine  i}i|. 
Let  us  assume 


The  boundary  condition  (2)  for  function  ,  gives  P^,  as  condition  to  S 

dF,  o 

^=2yi 

And  hence  the  problem  consists  of  determining  the  harmonic  function  P^  -  0), 

which  on  the  lateral  surface  of  the  cylinder  S  satisfies  the  condition 

^=0  (\) 


and  on  the  bases  of  the  cylinder  ^ 
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where  the  minus  and  the  plus  sl^ns  are  taken  from  the  lower  (U),  and  the  upper 
(o)  bases. 

In  the  oylindrioal  coordinates  the  Laplace  equation  takes  on  the  following 

I  1  ^  I  _L 

Hr*  '  r  dr  ^  r* 


Let  us  assume 


/’,  =  8in0  2C./?.shlA.(z- 


where  are  oonstants,  is  the  function  of  r  only,  is  the  ooordinate  of 
the  oylinder  upper  base  S,  is  the  ooordinate  of  the  lower  base.  Prom  the  La¬ 
place  equations  for  it  follows  that  the  functions  R^  must  satisfy  the  equa¬ 


tion 


dm. 


ar* 


+  T4r  +  (^:-?)"-=“ 


By  introduoing  new- variables  C  = ,  we  equate  this  equation  to  Bessel's 
equation 


;2  ^  C  rfC  \ 


d 


By  the  sense  of  the  problem  under  consideration  with  r  -  0,  one  should  get 
R  *  0,  hence  along  the  z  axis  of  the  cylindrical  cavity  S,  the  liquid  must  not 
be  of  infinite  values  for  the  velooity,  but  limited  with  r  ■  0,  the  magnitude 
, should  be  .  r~‘5P,/5Ci. 

Therefore,  the  fo^owing  should  be  accepted 

ft 

^  =  (Q  =  J  cos  (i  cos  f))  sin*  0  rfO 
0 

where  is  Bessel's  integral  of  the  first  order  and  first  degree.  The  bound¬ 

ary  oondition  (3)  on  the  lateral  surface  S  gives  the  relationship  dFJdr  =  0,  or 
dHjdr  =  0,  ,  or  better 


i/C 


=  0 


N.  £.  Zhukovskiy  determined  the  first  roots  of  this  equation 
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c,=  1.8412.  r>.43ir),  s.r)3i>3 

^4  =  11.7060,  !:s  =  14.8633,  =  18.0155 

The  values  of  the  constants  ,  are  deterained  by  the  formula  x  =. 


The  boundary  oondltlons  (4)  for  the  lower  and  upper  bases  gives,  the  fol¬ 


lowing  relationship 


2C>./,()..r)chX,/j  =  2r 


Since  ^  ■  2h,  Hence,  by  using  the  well  known  formulas, 

■  10,  CCJIH  H  ^  m 

I  (V)  M  r  dr=  _  1 . ,  ^ 

;  1— y-j — l-XilX.®)!*  ,  ecjiB  »  =  m 


the  constants  C  are  determined: 
n 


c.= 


_ f  J ,  (\,r)  r*dr 

(kna)P  c\^  iKh)  3 


Hence 


“  S  qctr,  ^  istsw''- 


Let  us  compute  ,  [ij*^  If  we  take  into  account  the  boundary  conditions 

(3)  and  (4),  we  will  have 

=  pf jj  (y*  +  **)  *  +  P  JJ 2y  do-Af  (z,-z,)  If  yVo 


where  C  is  the  base  of  cylindrical  cavity.  By  substituting  the  calculated  val¬ 
ue  of  the  function  we  obtain 

H,..  -““'[v-sfS  («) 


1IC^9/V 


where  ■A/  =  2irpa*/i  is  the  mass  of  the  liquid,  filling  the  inside  of  the  shell. 

This  formula  coincides  with  the  formula  of  N.  E.  Zhurovskiy  with  -  h. 

From  the  symmetry  of  the  cavity,  with  respect  to  the  a  axis,  there  follows 
•jhe  equality  =  •  It  is  also  difficult  to  determine  that 

iijl,,  ()i,]=0.  .  Let  us  designate 

Prom  this  the  kinetic  energy  of  the  rotating  motions  of  the  missile  that 
is  filled  with  liquid  will  he 

2T={A'-\-A*)(p*  +  q^)-^C'r^  . 

The'  stability  problem  of  rotating  motions  of  such  a  missile  is  parallel  to 
the  classic  problem  of  the  stability  of  a  simple  missile  with  the  moment  of  in¬ 
ertia  of  A  ■  B  -  A'  +  A*,  C  «•  C.  Hence  (3)  one  can  immediately  write  the  sta¬ 
bility  condition  of  rotating  motions  of  the  liquid  filled  missiles  on  flat  tra¬ 
jectories* 

CV-4(A'-H^*)pl>0 

For  an  ideal  incompressible  liquid,  and  for  the  same  type  of  missile  as 
studied  above  (C,  A'  ■  oonstants)  and  for  the  same  conditions  of  launching 
(r,  (A  -  oonstants)  the  density  of  liquid,  in  increasing  the  teisn  A*,  has  a 
negative  effect  upon  the  conditions  of  stability. 

3.  Cylindrical  Cavity  with  One  Flat  Diaphram.  In  order  to  increase  the  axial 
moment  of  inertia  of  a  missile  filled  with  liquid  and  thus  increase  the  stabil¬ 
ity  of  rotating  motions  of  the  missile  on  its  trajectory,  the  diaphrams  are 
placed  in  the  cylindrical  cavity. 

Let  us  study  the  influence  of  one  flat  diaphram  set  radially  in  the  cavity. 

A  missile  with  one  diaphram  does  not  possess  axial  symmetry.  Hence,  to 
simplify  the  oaloulatlons  for  the  x,  y,  z  axes,  we  will  ohoose  axes  tied  with 
rigid  shell,  with  the  origin  in  the  center  of  gravity  0  of  the  filled  missile 

so  that  the  plane  y  ■  0  would  be  the  plane  of  the  diaphram. 
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Projections  of  instanti  neous  absolute  angular  velocity  of  the  missile  on 
this  axis  will  be 

p  —  ^'cosMi  j  a' sin 0) cos ^ 
q  =  —  sin  ui  -j-  a"  cos  o)  cos  ^ 
r  =  ut'  —  a'siiipl 

Kinetic  energy  of  rotating  motions  of  the  rigid  shell  will  be 

2T'=A'p^-{-n'q^-{-C'r^ 

where  the  constants  A',  B',  C  represents  the  moments  of  inertia  of  the  missile 
with  respect  to  the  axes  x,  y,  z.  The  kinetic  energy  of  the  liquid  will  be  ex¬ 
pressed  by  the  general  formula 

2r*=  2  t'l'o 'W 

It  is  necessary  to  determine  the  function  <]>(.  .  The  function  'j's  was  found 
by  Stock. 

( 

We  will  study  that 'part  of  the  cavity  which  is  divided  by  the  diaphram  and 
in  which  (  (i  =  rcos9,  y  =  rsinO)  )  the  angle  0  .varies  from  0  to.it.  .  Let  us 
assume 

The  condition  (2)  on  the  boundary  of  the  cavity  gives  the  relationship 

dFJdn  =  —  2ya. 

Prom  this,  on  the  wall  of  the  cylinder  (  (0<9<k))  we  have 

dFJdr=  —  o  sin  20  {r  =  a) 

on  the  diaphram  (9  =  0,  9  =  it) 

dFJdr^  =  0  (8) 

at  the  sections  of  the  lower  and  upper  bases 

dFJdz  =  0  (9) 

We  will  try  to  find  function  in  the  form 

CnF'cosnB 

( _ _ _ _  n 
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The  boundary  condition  (9)  is  satisfied  immediately,  sinoe  P  does  not  de¬ 
pend  on  Z  condition  (8)  is  also  fulfilled,  since  the  walls  of  the  diaphram 

1 

y  >  0  we  have  9  ■  0,  =  •  If  we  use  the  formula 


For.  the  half  (  y^O)  of  the  circular  oavity  under  consideration  we  oan 


calculate 


14’r  'I'jl  — «*  S  n(n*  — 4)*)  ("  =  1.3,5,...) 


where  as  before  M  =  2«^h.  *  show  that 

(tj.  W— 0.  W  = 

Let  us  calculate  the  function  'h- .  Let  —  — y*:  the  boundary  con¬ 

ditions  (2)  give  us  the  following  relationship  on  the  boundary  of  the  oavity 
under  consideration 

d/’,/dn  =  2y7 

In  other  words,  on  the  lateral  surface,  sinoe  7  =  0,'  it  should  be 

dP,/d/i=0 

and  on  the  bases 

dFJdn=T2y  (^2) 

where  the  plus  and  minus  signs  are  taken  for  the  lower  and  upper  bases,  respec¬ 
tively. 

Inside  the  oavity  the  function  P^  must  be  harmonic  ^  Pj^  -  0}  in  cylindri¬ 
cal  coordinates  it  must  satisfy  equation  (^) 

Let  us  assume 
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where  R 

nun 

tion 


denotes  only  the  function  of  r.  Function  R  must  satisfy  the  equa- 

nm 


If  we  introduce  a  new  variable  C  =  \*,/‘,  then  it  will  be  the  familiar  Bes- 

ael  equation  ,  f. 

“5C*“  +  r“^SC  cw  - 

For  this  problem  we  must  take  « 

=  /,  (0=  {  cos  (C  sin  z  —  m)  dz 

0 

where  /.(C)  denotes  the  Bessel  Integral  of  the  nth  order  and  first  degree. 

On  the  surface  of  the  cylinder  (  (r  =  a,  0<0<ii)  )  from  boundary  condition 

(11)  we  derive  the  equations 

d/,(C)/«iC  =  0  (c  =  x,«ii) 

■  I 

determining  the  eigen  values  of  >nm-  The  boundary  condition  (ll)  on  the  dia- 
phrara  (y  -  0,  6  -  0,  6  ■  Tf  )  is  immediately  satisfied.  In  accordance  with  the 

condition  (l2)  for  the  upper  and  lovrer  bases  we  have  (Z^  -  «  2h) 


2  cos  ch  (Kmfi)  =  2r  sin  0 


(13) 


For  the  value  0  within  the  interval  (O,  )  the  function  sin  d  may  be  ex¬ 

panded  into  Fourier  series  in  accordance  with  the  cosines. 

By  substituting  this  series  in  formula  (l3)  and  equating  the  coefficients 
of  Cos  $  f  VQ  will  obtain  for  n  ■  0,  2,  4»  *  •  , 


y.  (>-,m  r)  ch  (Kmh)  =  ”  7?^ 

m 

and  in  the  first  part  (of  equation)  zeros  with  n  -  1,  3»  5»  ••• 


(14) 
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j  rJ^  U^kT)  y,  (/,„r)  dr  = 


Kla-i 


2  [■/»  (*ii»i<*)I‘  npH  ^nzxXat 


the  following  terras  of  the  oonstanta(are  determined)  (n  -  0,  2,  4)  •••) 


C. 


16 


mm  Q 


7J  ry,Q,„r)dr 


For  odd  terms  n,  the  constants  C  are  zeros. 

'  nra 

Substituting  these  constant  terms  in  the  function  and  the  latter  in  the 
already 'used  expression  (2)  / [<}>,,  »  after  calculation  we  obtain 

22L  V  >■?!>■  _ (  V  -  -  f  rV.  (K»r)  d/-)*] 

+  K*a*A  l)(n*— D*  \ ^  '  ■* 

Now  we  will  find  function  .  .jjj.  With  this  in  mind  let 

■|'2  =  «+^2 

The  boundary  oondi-y.on  on  the  surface  of  the  examined  oavity  will  be 

dt\!dn  —  —  2.T7 

This  means,  that  on  the  lateral  surface  of  the  oavity  I  (-[=0)  )  we  must  have 

dF^jdn  =  0 


and  on  the  bases 


:  dFJdn  =  +  2x 


where  the  plus  and  minus  signs  are  taken  for  the  lower  and  upper  bases  respeo- 

I 

tively. 

We  will  seek  function  P^  in  the  form  of  the  series 
/■,  =  2  C.  cos  0/,  (Xy)  8h  X.  [z  - 

It  is  on  that  part  of  the  surfaoe  of  the  cylinder  which  enters  the  bound¬ 
aries  of  the  half  of  the  oavity  under  consideration  virtue  of 

boundary  condit  ions  we  have  the  equation 
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I  _ Q  ■ 

da 


determining  the  constant  terms  • 

On  the  flat  diaphram  the  boundary  condition  is  satisfied,  since  there 
sin  ^  ■  0. 

Prom  ‘boundary  condition  of  the  bases  of  the  oavity  we  get 

^C,J,{\,r)\,ch{\.h)  =  -2r 


Hence 

Therefore 


C.= 


4a* 


•  ~  C,  (U*  —  1)  ch  (K>>)  J\  (C) 


cos  G  2  ^  _  j- j  j  ^  j  sh  (C. 


2i  — Jo- 


Substituting  this  expression  for  in(  4'*-  4'i  )  after  calculations  we  get 

;  A/  r  Jo*  -t-Jo*..-f*«*  I  5  a  V _ 1 _ th  (r  A)] 

Itj.  W  — tL - 3  ^40  '‘a  A  V’  flJJ 


For  the  oavity  under  discussion  (  "'I'l*  t*  )  -  0, 

The  moments  of  inertia  of  the  accompanying  mass  of  the  liq.uid  filling  the 
cylindrical  oavity,  which  is  divided  by  one  flat  diaphram,  will  be 
^*=2[t„  til.  B*=2[tr  til.  C*  =  2[t3,  tsl 

The  kinetic  energy  of  the  rotating  motions  of  the  missile  together  with 

I 

the  liquid  filling  it,  will  be 

2T  -  Ap2  -  Bq2  -  0^2 

where 

A  -  A*  +  A*,  B  «  V  +  B*,  C  -  C  +  C* 

From  this  we  know  that  the  ellipsoid  of  the  inertia  of  the  missile  and  accom¬ 
panying  masses  will  be  triaxial  in  this  case.  Such  a  case  of  an  elongated  mis¬ 
sile  with  solid  filling  has  not  been  completely  investigated  for  its  stability. 

It  is  necessary  to  examine  it  separately. 

I 
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With  this  in  mind,  and  using  for  the  Lagrange  variables  the  angles 

we  set  up  equations  of  the  rotating  motions  of  the  missile  in  the  form 
of  the  Lagrange  equations*  In  doing  so,  we  take  as  the  generalized  foroes,  as 
was  explained  in  (l) 

^  ^.  =  |i8inacosp,  =  jjl  cos  a  sin  ^„  =  0 

The  obvious  form  of  Lagran^  equation  in  this  oase  will  be 


C  ^  (u)'  —  a'  sin  P)  =  (yl  —  B)  ( —  p'*  -j-  a'’cos‘'‘  p)  cos  «>  sin  lu  -|- 
4  {A  —  //)  a'p' (cos’^  to  —  sill*  to)  cos  p 

(p' cos  to  4'='"  sin  to  cos  p)  sin  u)  cos  p  •- 

4-  B  ( —  p'  sin  to  4-'  o!  cos  u)  cos  P)  cos  oi  cos  p  — 

—  C  (to'  —  a!  sin  p)  sin  p)  =  (a  sin  a  cos  p 

^  (/I  (p'  COSO)  4- a' sin  to  cos  p)  cos  to  —  B  ( —  p'siii  to  4-  a'cosoicosP)  sin  toj  = 
=  —  [  •^  (P'  cos  to  4-  sin  to  cos  p)  a'  sin  oi  sin  p  4“ 

4-  B  ( —  p'  sin  to  4-  2'  cos  to  cos  P)  a!  cos  to  sin  p  4“ 

,  4-  C  (to'  —  a!  sin  p)  a'  cos  p]  4“  cos  a  sin  p 


These  equations  have  a  particular  solution 

I 

0)'  =  to,',  a  =  0,  .  p  =  0 

Let  us  take  this  particular  solution  for  the  undisturbed  motion  and  try  to 
investigate  its  stability  in  the  first  approximation.  Let  u),  =  to,'<,  to  — 0),  =  ;. 

The  equations  in  the  variations  for  this  vmdisturbed  motion  are  of  the  form 

(7r=o  ' 

A  (p"  cos  t..,  4-  a'  si n  0.,)  4-  (/I  —  fi)  <  (—  P'  sin  to,  4-  a  cos  to,)  = 

=  —  Co),'  ( —  p'  sin  0),  4*  a'  cos  to,)  4- (P  cos  to,  4-  a  sin  to,) 

B  (—  p" sin  u),  4-  a"cos to,)  4-  (/I  —  5)  to,'  (p' cos  to,  4-  a'  sin  to,)  = 

=  Cto,'  (P'  cos  U),  4-  a'  sin  to,)  4- 1*  (—  p  sin  to,  4-  a  cos  to,) 


MCI/^9/V. 


-13= 


This  system  presents  equations  with  the  periodic  ooeffioients  (*4} 

As  Lyapunov  established,  suoh  systems  of  differential  equations  can  be  trans* 
formed  into  a  system  with  constant  coefficients,  without  changing  the  problem 
of  stability.  The  first  of  the  equations  is  of  no  particular  interest;  it  ex¬ 
presses  the  obvious  properties  of  the  angle  hP  when  the  given  terms  are  affected 
in  the  first  approximation.  In  order  to  transform  the  last  two  equations,  we 
will  introduce  new  variables 

u=asinu)o4-?coso.o.  i»  =  aco8u),— pisinu.. 

The  determinant  of  the  transformation  is  different  from  zero  and  ■  -  1. 
Hence. 

a'sinu)g-|-p'co8ujj  =  u'  —  to^'v  ■ 
a.'  cos  (Bj  —  sin  u)j  =  w'  -|-  uj^'u 

a"  si  n  (Bq  -j"  P”  =  u"  —  2u>^v'  —  u),'*u 

a"  cos  u)o  —  P"  sin  u>j  =  v"  -j-  2v>^'u'  —  lu^'V 

It  follows  then  that  the  equations  in  the  variations  with  new  variables 

will  be  >■  i4(u"  — 2(1), V  — «<>(,'*«)  =  + 

B  (v"  -f-  2to,'u'  —  «>o***’)  =  "i*  ^“’o’  — ‘"o’*’) 

or  better  + [(C « +(C' — e)u,,V  =  0 

Bv" 4- [(C  —  A) .B,'*  —  ^\v  —  {C  —  A  —  B)  =  0 

The  characteristic  equation  of  this  system  of  differential  equations  with 
constant  coefficients  is 

+  { [  ^5  +  (A  —  C)  (fi  —  C))  «B,'*  - 1*  (i4  —  fi)}  X* + 

The  conditions,  that  the  roots  of  this  equation  are  purely  imaginary,  are 

expressed  by  the  following  three  inequali tries"  . 

I(C _ (i)  [(C - 5) iB,'*  - ,i) > 0 

I(A5  +  (^-C)(5-C))(b;*-ix  (^  +  5))*-  I 


We  note,  in  the  last  condition  when  A  -  B,  leads  us  to  the  well-known 
Mayevskiy  theorem  of  the  stability  of  the  rotating  motions  of  a  missile;  in 
this  case  the  first  two  inequalities  will  be  satisfied. 

4.  Cylindrical  Cavity  vd.th  Two  Crossed  Biaphrams.  Let  us  examine  the  case 
where  two  dlaphrams  are  orthogonally  placed  with  respect  to  each  other  in  the 
cylindrical  cavity.  Since  the  ellipsoid  of  Inertia  of  the  missile's  rigid 
shell  will  be  the  ellipsoid  of  rotation  around  the  axis  Z,  the  axes  of  the  co¬ 
ordinates  bound  to  the  dlaphrams  are  so  chosen,  that  their  equations  will  be 
X  ■  0,  y  -  0  and  later  in  the  dynamic  part  of  this  paragraph  we  will  change 
to  the  axes  used  in  (2). 

Function  tjij  was  determined  by  Stokes  taking  the  limit  in  calculating  the 
value  of  one  unknown.  There  is  no  point  in  changing  Stokes'  calculations,  and 
we  will  accept  the  expression  (  .'j's.  ’I'j,  )  as  , derived  by  Zhukovskiy.  For  the 
cavity  under  discussion 

ll'i.  I'i'r  ’l'sl==0.  (I'l.  1*11  =  0 

and  it  is  evident  that  after  finding  function  the  problem  of  finding  func¬ 
tion  <{>f  does  not  arise,  since 

•  (1*1.  Ill  =  11*1.  I*,) 

For  the  purpose  of  computing  <!»,,  ,  let  I*,  =  — 


Boundary  conditions  (2)  give 


On  the  lateral  surface  of  the  cavity  0^6^  l/2ir, ,  separated  by  the  crossed 
diaphrams  Y-  and  it  should  be 

dn  ® 


and  on  the  bases 


where  the  minus  and  plus  signs  are  taken  for  the  lower  and  upper  bases. 
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We  shall  look  for  the  hannonio  function  in  the  form  of  series  expan¬ 
sion  (n  -  0,  2f 

f^  =  y  C,„ cos  (>^„r)  sh  ( Km —  i --)) 

N,  m 

On  the  surface  of  the  cylinder  limiting  the  cavity  0^0^'/, «,•  under  dis 

oussion,  we  have  in  accordance  with  boundary  conditions 

O  kihO) _  „ 

da 

These  equations  determine  the  value  r.  The  boundary  conditions  on  the 

I 

diaphram  will  be  satisfied,  since  when  n  is  even  on  the  diaphram  we  will  have 

sin‘n^  «,0-.  The  boundary  conditibnsi  on  the  upper  and  lower  boundaries  bases 

give  ... 

•  2  COS  «0/.  ()..„r)  cli  (/.,„//)  —  2r  sin  << 


By  takingO<  0  ^‘/j'used  in  series  expansion  of  sin  ^  in  the  preceding  para 
graph  we  obtain  the  established  values  of  Cmm  when  (n  ■  0,  2,  4, •••)• 

r*. 

64p  v  1  X,.„  111  (>,„  A)  /  1  i’  9/  ,,  -)ffA 

From  this  the  moments  of  inertia  of  the  accompanying  masses  of  the  liqu^ 
that  fill  up  the  four  sections  of  the  cavity  separated  by  the  orthogonal  dia- 
phrams  will  be 

A*  -  B*  -  (•;„  ■i.i.  c*  -  4  [  ['>3.  ’fsl’  ]. 


The  ellipsoid  of  the  inertia  of  a  missile  that  is  filled  with  an  ideal 
liquid  will  be  the  ellipse  of  the  missile's  rotation  along  its  axis  of  propa¬ 
gation.  The  problem  of  the  stability  of  rotating  motions  of  such  a  missile 
coincides  with  the  classical  problem,  and  hence  we  can  immediately  write  the 
condition  for  stability 

+  — 4(4'+4*)  ii>0 
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On  the  left  side  is  the  second  degree  polynomial  of  the  liquid's  density 
,  Hence  the  stability  of  the  cavity  will  be  determined  by  the  disposition 
of  the  real  roots  of  this  polynomial. 
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